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Abstract 

We derive coupled propagation equations for ultrashort pulses in a degenerate three-wave mix- 
ing process in quadratic media, using approximations consistent with the slowly evolving wave 



0,| approximation [T. Brabec and F. Krausz, Phys. Rev. Lett. 78, 3282 (1997)]. From these we de- 

o ■ 

rive an approximate single-field equation for the fundamental field. This document expands upon 



mathematics used for work submitted by the same authors to Physical Review Letters. 



I. INTRODUCTION 

The following derivation has two main parts. First, Maxwell's equations are reduced to 
a set of coupled propagation equations for degenerate three-wave mixing in quadratic media 
with less restriction on the minimum pulse duration than required by the Slowly Varying 
Envelope Approximation (SVEA). To do this we use a set of approximations consistent with 
those used by Brabec and Krausz in their derivation of the Nonlinear Evolution Equation 
(NEE) |l|. These are known as the "Slowly Evolving Wave Approximation (SEWA)". Sec- 
ond, we collapse the coupled propagation equations to an approximate single-field equation 
for the fundamental field using perturbation methods. 

II. DERIVATION PART lA: GENERALIZED EQUATIONS 

We begin with Maxwell's equations for the electric field vector E in the Fourier domain, 
including linear P*^^) and quadratic nonlinear P*^^^ contributions to the electric polarization, 

{£, + Vi) E(r, .;) + ^E(r, u) = -^[P« (r, u) + p(^)(r, .)]. (1) 

Next, we expand the electric field and polarization into waves propagating at the funda- 
mental (FF) and second harmonic (SH) frequencies. In performing this step we assume the 
FF and SH spectra do not overlap. Here, ei and 62 denote the polarization of the FF and 
SH fields, respectively. 

E(r, uj) = Etoi(r, uj) = Ei{r, uj)ei + E2(r, uj)e2, (2) 

P«(r,^) = x?{^) ■ Etot(r,a;) + x? (u;) ■ Etoi(r,a;), (3) 

P(2)(r,^) = Etot{r,u;) ■ xf\u;) ■ Etot{r,u;) + Etot{r,u;) ■ x?\uj) ■ ^tot{r,u;). (4) 

The subscript j in Xj i^) denotes the k^'^-ordei susceptibility that gives rise to a FF or SH 
wave (for j = 1,2, respectively). 

Defining the linear indices of refraction [noj(c(j)] = 1 + 4:TTXj (uj) and wavevectors 
[^^(co')] = [r7,o,j(^)] itP'jc^ in the usual fashion and rearranging (1), we obtain a set of 
coupled-field equations, 

a2 \ /I 2 

+ Vl ) ^i(r,cu) + [A;i(cu)]'Ei(r,cu) = -^Etoi(r,cu) • ^^\uj) ■ Etotir,uj), (5) 



Q^2 ' -X ,^.V^>-/ ' L-XV-/J ^.V^,^y ^2 



^ + Vl ) ^2(r,cu) + [k2{uj)f E2{r,uo) = -^Etoi(r,cu) ■ xf M ■ ^tot{r,uj). (6) 



We now introduce the complex envelopes Aj and their Fourier transforms, 

Ei(r,t) = Ai(rx,z,t)e^(^^(^o)""'^°*) +C.C. ^ — > Ai(r^, z,cu - cJo)e*('^i('^«)") + c.c, (7) 

E2(r, t) = A2(rx, z, t)e^('=2(2^o)^"2^ot) + ^^p. ^ > A2(rx, 2, cj - 2cJo)e*(^^(^^"^") + c.c. (8) 



Substituting envelopes (7), (8) into equations (5), (6) and defining A A; = 2ki — k2, we 
obtain 

2iki{uJo)— + — - [ki{uJo)f + V\] Ai{r±,z,uj-iJo) 

+ [ki{uj)f Ai(r_L, z,uj- ujq) 
^'''''-X?l-.„Al{r^, z,oo + uo)A2ir^, z,u - 2c.o)e-^^^ (9) 



2 'Mu>o-(^o='^0 



^2 



2ik2{2uQ)— + -g^- [k2{2uo)f + Vi 1 A2{r^, z,uj - 2cc;o) 

+ [k2iuj)f A2{r±, z,u- 2uJo) 
= -^x£V.o=2.„^?(r±, z,u;- u;o)e^^'\ (10) 

(and complex conjugates). Henceforth, we'll ignore the conjugate equations. 
As usual, we expand wavevectors kj{ujj) about ujj {uji = ujq, 002 = 2ujq): 

kjiuj) = kjiujj) + k'jiuj — ujj) + Djijjj — LJj), (11) 

with k] = ^U^ and Dj{uj - uJj) = Y.n^^ '^{uj - ^^O"- 

We will now assume the quadratic nonlinearity is instantaneous, i.e. non- instantaneous 



effects are negligible 






) , and perform the reverse transform: 



«*4-^-*?^^i 



+ 



ki + tk[— + D,{t) 



1 2 



57r 
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-^ [ujo + i—] Al{r±,z,t)A2{r±,z,t)e 



'iAkz 






) 



(12) 



+ 



k2 + ik^^ + D2{t) 



A2{rj_,z,t) 
A2{rj_,z,t) 



-^(2c.o + ^|) A?(r„.,t)e^^- 



(13) 



where Dj{t) = X]^2 "^ (^^)"- Note: o;^ = (cu — a;o)^ + 2(a; — ljq)ljo + cjg, and thus 



UJ 

^-'[u'GiLu)] = -^^ + 2iLUo^ + ui = (luo + t§ifGit). Also, u^ = (cu - 2iUoy + 



d^Gjt) 



2{Lj-2LOo){2uJo) + {2uJoy. 

Next, we transform coordinate systems to co- moving variables in the frame of the FF 
wave, T = t — k[z, z' = z. Partial derivatives become 



d_ _ d_ 



d d , d 

dz dz' '^ dr ' 



> d' , ,j2 d' 



(14) 
(15) 
(16) 
(17) 



dz^ dz'^ dz'r dr^ 

(We drop the prime on z for notational simplicity) and the equation for the FF becomes 



2*i^-*;^VfS-2*;^^ ■"'' 



dz dr 

+ 



dz"^ 



dzdr 



+ k 



12 



kl + V] 



kf + 2ikik[^ + 2kiDi - A;f ^ + 2ik[D^^ + Dl 
or OT^ or 



Ai 

Ai 



Stt 



c^o + i 



d_ 



X2aJo-'^o'^l'^2e 



(18) 



Combining terms, the equation for the FF simphfies to 

2iki (l + i 



.k[ d\dA, _^2^^_d'A, 



ki dr J dz 



The equation for the SH becomes 



a.= -2*-«-l' + 'f|l^'-^'^- 



^7r^oxS„-..„ 



1 



i d 
ojq dr 



AlA^e 



-i!S.kz 



(19) 



liki 



dz ^ dr 



+ 



n2 a2 a2 

^ 2k' — h k'^ — 

dz^ ^dzdr ^ dr^ 



+ 



kl + 2ik2k' — + 2k2D2 - H^tt^ + 2ik'^D2^ 



kl + vl 
d 



■dr 



di 

47r / .d 

■^ 2-0 + ^^ 



dr 



Di 



A2 
A2 



p) 



Aa;o+i^O 1 



2AAkz 



(20) 



Since our coordinate system is in the frame of the FF wave, fewer terms will cancel when 
we combine terms in the SH equation: 



2iko 



k[ d\ dA2 „2 . d^^2 



/c2 dr J dz 



-VU, 



2koDo 1 



dz^ 
-2ik2{k'^ - k[ 



■^2 d\ , ^2, 

,aA2 ,,,2 ,,2^5'^2 



dr 



1 + 



(^1 ^2 ) 



dr^ 



d 



a2 lAkz 



(21) 



c^ \ 2uJo dr ^ 

The fifth term on the RHS is the familiar group velocity mismatch (GVM) term, while the 
sixth term on the RHS is unfamiliar. Also notable, unlike in the FF equation, the prefactor 

dA 
dz 



to the §7 term has mixed wavevectors by field ( ^ ) • Also, the GVM term has no prefactor 



r-^ 2 



Assuming we'll want to use the approximation Vp c:^ Vg ( ^ ~ -^ for FF, and -^ _ j, 
for SHJ to simplify the equations in accordance with Brabec and Krausz's SEWA, we use 
some algebra to rearrange terms in the SH equation: 



^1 ^2 ~ (^1 ^2/(^1 + ^2) 

= \k^ — k2){ki — A;2 + 2^2) 



(K 



k'2? 



2^2(^1 ^2. 



(22) 



Thus, 



Also, 



2ik2(k'2 — k[ 



,aA2 



dr 



{k'i 



12 



"2 ) Q^2 



d^A, 






-2ik2{k'2 - k[] 



^ ^ ^k[ d\ dA2 
k2 dr J dz 



1 + — — 



dA 



(^2 ^1) 



d'A2 



dr^ • 



Xk'^-k'^ + kj) d \ dM 
k2 dr J dz 

k'2 d .{k'2-k'i) d\ dA2 



l+iT^^ 



k2 dr k2 

dA 



dr J dz 






ko dr J dz 



2 ^{k'2-k[)d^A2 



dzdr 



(23) 



(24) 



Using (23) and (24), the SH equation becomes 
.k'o d \ dA2 



2ik2 \ 1 + i 



ko dr / dz 



-V^As- 



d^A2 ^,,, ,,,92^2 



dz'^ 
2ik2{k2 — k[) 

2k2D2 (l - 



2(^2 ki^ 



dzdr 
k' d\ dA2 



(^2 ^1) 



d^A2 

dr^ 



k2 dr J dr 



;^— 1 A 

kodr 



1 + 



2 - ^^2^2 



d 



A{e 



2JAkz 



(25) 



c^ \ 2ljq dr ^ 

where the unfamihar third and fourth terms on the RHS are due to our choice of co-moving 
variables in the frame of the FF wave. 

Dividing through by the prefactor to the ^ term in each equation, our two generahzed 
propagation equations for FF and SH waves with instantaneous quadratic nonlinearity are 
(with 5 = k[-k'2) 

dAi 1 / k\ dY'2 , 



2.,(2) 



4vr^oXa 



1 + 



i d 
LOQ dr 



1 + i 



kic^ 

_!_ / .k{d_ 

'2k[ y^%d^ 



K d 



fcl dr 



■AlA2e 



-iAkz 



dz'^ 



+ Dl 



^1, 



(26) 



.dA, 



+ 



1 



2ko 



dz 2h 



kodr 



■ k'2 d 
kodr 



^^^0Xt^0+a;o 



VIA2 - i6 



1 + -J-^ 



dr 

2 



DoA 



2^2 



k2C^ 



-^ ^ ^ fc2 (9r 



a2 iAkz 



ffi Q2 Q2 



dz"^ 



dzdr 



dr^ 



A, 



(27) 



III. DERIVATION PART IB: THE SLOWLY EVOLVING WAVE APPROXIMA- 
TION 



Here we make the first approximations (other than that of instantaneous quadratic nonlin- 
earity) that place a restriction on the minimum- duration pulses our equations will accurately 
model. Brabec and Krausz point out that the conditions of the SVEA, 

dAi 



dz 
dA, 



dr 



<^ kj\Aj\, 



<^ujj\Aj\, 



(28) 
(29) 



can be relaxed by replacing (29) with 



kj 



Ujkr 



k, 



«1, 



(30) 



(i.e. phase and group velocities are approximately equal), resulting in the NEE. Brabec and 
Krausz called this combination of approximations, (28), (29), the SEWA, since it can be 
written compactly as 

-g^ «k,\E,\. (31) 

Equation (31) illuminates a key difference between the SVEA and SEWA: the SEWA requires 
the field rather than envelope to not change significantly as it propagates over the length 
scale of a wavelength, and the approximation no longer directly restricts the bandwidth of 
the wave to be smaller than its carrier frequency (or the pulse duration to be much longer 
than a single optical cycle) (29). 



Condition (30) will be satisfied for our coupled FF and SH waves if both 



Ulo 



fcl 



and 



\2c<;o 



fc2 



are true. Applying these conditions to equations (26) and (27), the equations 



reduce to 






fclC^ 



i d 



AlA^e 



lAkz 



(32) 



.dA2 1 
^ a^ 2A;2 



1 
2A^ 



9 



2ujr, dr 



VIA2 



dA 
i6-A + D2A2 



^TT^O^XS+^O 



fcsc^ 



t2 

i d 



1 + 



9r 
9 



2ujq dr 



2iAkz 



Afe 



-1 



26- 



d^ 



dzdr 



+ 6' 



d 



i2 H 



9r2 



A, 



(33) 



92 A, 



"aP" 



< 



aA, 



2ct;o dr ^ 

where the -q^ terms are neghgible by condition (28), since it also imphes ^ 
We also have eliminated the D'^Aj terms, assuming all dispersion terms of order -^ and 
higher will be negligible. (* This is a necessary assumption of NEE, NLSE and all first-order 
envelope equation derivations, and we similarly use it here.) 

The equation for the FF closely resembles a ID temporal version of Brabec and Krausz's 
NEE, with quadratic nonlinear term replacing cubic. However, the SH equation contains 
the two new terms 



S d^ 



and 



52 d^ 



k2 dzdr 2k2 dr"^ 

Our next task is to consider if either of these terms must be eliminated on the basis of being 
self-consistent with the approximations we've already made. 

Using a similar approach as P. Kinsler uses in his derivation of the GFEA equation jj] , we 
find it is self-consistent at this point to eliminate the 



i-TTir term, as demonstrated below. 

fc2 OZOT ' 



We can write each of our coupled equations in the form 

.dA, 



dz 



+ [LHSj] = [RHSj]. 



For this equation, the condition 



is true if and only if 



\RHSj\ < 



\LHSj\ - 



dA, 



dz 
dA, 



dz 



(34) 



(35) 



(36) 



is true as well. (Note, this simply amounts to \LHS\ ^ \RHS\.) We have already made the 
approximation (28): 



dA, 



dz 



^ ^jl^il- 



Combining (28) and (36), the following condition must hold as well 



\LHSj\ < kj\Aj\. 



(37) 



We have already eliminated -q^ terms by adding a ^ to each side of (28). Thus, we 
have assumed (following (37)), 



d , 

^\LHSA < 
oz 



^ dz 



Since 



it follows that 



.8 A, 



dr 

S d^A2 



< \LHS2l 
dA2 



(38) 



(39) 



ko dzdr 



< 



dz 



(40) 



proving that the elimination of the -g^ terms is consistent with the elimination of -q^- 
However, using the same type of argument we cannot determine that 



6^ d'A2 



ko dr"^ 



< 



dA, 



dz 



(41) 



5292^2 



and it is not consistent to eliminate the t-^k^ term. 

fe2 OT 



IV. THE NONLINEAR EVOLUTION EQUATIONS IN QUADRATIC MEDIA 



Thus, after using the SEWA our propagation equations become 

■9Ai 1 / i dY^^o , 



kic^ 



i d 

ujq dr 



AlA^e 



iAkz 



0, 



(42) 



.dA2 1 / id 

dz 2k2 \ 2(jJo dr 



.2 . , .,d^A2 



VU2 + S' 



dr^ 



- ^6^-^ + D2A2 + ^^"o^-U 



dr 



k2C^ 



1 + 



d 



2uJo dr 



2iAkz 



Aie 



(43) 



which are the Nonhnear Evolution Equations for coupled FF and SH in a quadratic medium. 
After transformation to a nondimensional coordinate system, ^ = {5/tq)z = z/Lgvm, 
s = t/tq, where tq is the initial temporal half- width, the NEEs become 



^^ + ^ 1 + IT ^i«i + ^i«i + 1 + TT «I«2e-^^^ = 0, (44) 






i^ + -7r 1 + 



9^ 2 \ 2u;oTo 9s 



P2 c^s' 



2 



da2 ^ 

- «^- + "^202 

as 



^ 1 1 ^ 2:^1 1 "?»■" - °- («) 



with definitions, 



dihhy/^c^ She 

-^i CiCli, C\ — - J—: J—: — , C2 



(3 = LcvM^k = ^^, Pi = — = —-, iy= J- = — — , 

ki oki TQk2 ro/c2 

V. DERIVATION PART 2: THE SINGLE-FIELD EQUATION 

n 

Using the same perturbation- method technique as Menyuk et. al. 3], we collapse coupled 
propagation equations (44), (45) to an approximate single- field propagation equation for the 
fundamental field (FF). 

We employ the transformations hi = ai/lP]^^"^ and 02 = 026"*'^^, and restrict our analysis 
of transverse dimensions to the temporal only. Furthermore, since dispersion terms above 
second order in (44), (45) will not give rise to significant terms in our derived single- field 
equation, we will ignore them here for simplicity. Doing this, we obtain 

,80,2 .80,2 0:29^0.2 z/ / i d\~ (9^a2 

' sf " ' ' ~ ' sF ~ y &^ ^ 2 V 2^9jJ "aJ^ 



10 



where aj = Lgvm/LgvDj = k'- /tqS. We can treat the fifth term of equation (48) as a 
perturbation to GVD and higher orders of dispersion, since they are of the same order in 
our expansion. Thus, we define ttg = '^2 — ^, and equation (48) becomes 



.80,2 n^ M2 a'^d^tti , |.| , _, , 



0. 



(49) 



d^ '" ^ ds 2 ds'^ ' I'-i ^ ' 2ujQTods^ 

In order to find an approximate expression for 02 in terms of Oi, we now expand the SH 
field in a power series in 1/|/9|, 



02 



00 

E0'2,n 
n=0 "^' 



and substitute this in the coupled equations (47), (49): 



dai tti d'^&i 

«-TT Z TT^ + 1 



d^ 



2 (9s2 



d 



uqTo ds 



0-20 + 



^21 , ^22 



+ 



1/31 m' 



(50) 



(51) 



.d 



a2i a22 , 

020 + 77^ + TTTTTT + 



1/91 I/9P 



/9 



020 



021 ^22 

WW 



— « 



+ 



d_ 
ds 



021 022 , 

020 + TT^ + TT^ + 



1 + 



i d 



2 9s2 



021 , O22 , 
O2O + TTT7 + TT^ + 



a?) 



2ujqTq ds 
Keeping terms to order 1/|/3| only, we get 



0. 



.dai aid'^ai 



di 



2 ds^ 



1 + 



d 



uoTq ds 



{a^a2Q) 



m w 



d 



uoTq ds 



0'lO'21 



(52) 



(53) 



da2o i da2i ^^ P ^ P ^ da2o i da^ 



21 



a'n d (120 cn'o d 0,21 



+ 1/91 1 + 



d 



0. 



(54) 



2 ds^ 2\p\ ds"^ ' ""' \' ' 2iUoTods^ 

Now we collect terms of like order in (53), (54) and use substitution to eliminate 020 from 
(54). From terms of lowest order (order |/5|), we obtain 

' ^^-- (55) 



P \ 2uJoTods 



At the next order (order 1), we have 

,dai aid'^ai 



dC 



2 ds^ 



d 



uJoTo ds 



(01020) = 0, 



(56) 



11 



^^^ - — «2i - ^^^ - — ^j::^ = 0- (57) 



.5a20 _ /5 -- _ •'9a20 _ 0:2(9^020 

Substituting (55) into (56), the order-1 equation for the FF becomes 



i |/9| „ <9|aip 1 \I3\ (daid\ai\'^ g^^^i 



Substituting (55) into (57), the order-1 equation for the SH becomes 

„.^ dai 1 fdaidai ^ d'^ai\ ^ „_ dai 
2«ai— — — — - + ai— — - - a2i - 2iai— - 

•yC 1^0 To V '^^ "^S OSOt, J OS 



^ 1 2cuoro 2 J 9s2 4cuoro Ss^ °- ^^^^ 



We can now ehminate %- terms by substituting (58) into (59), and the resulting expres- 
sion can be solved for 021. Finally, with both 020 and 021 known, we can write 62 to order 

1/1/31, 

02 ~ 020 + T-oT- (60) 

The resulting expression, containing many terms, is then plugged into our original equation 
for evolution of the FF (equation (47)), leaving us with our desired result: an approximate 
single-field expression for the FF. 

We now return to our original definition of the FF, ai, through the transformation ai — > 
ai/\(3\'^^/'^\ Our single-field equation now contains terms of order 1, 1/(3, l//5^, and l//9^. 
Of the more than 60 distinct terms in this NLSE-like equation, only a few higher-order 
correction terms are significant. We find, to order l//5^, 
dai «! d'^ai 1,12 1 1 12 9ai 



13 ujqTq \ ds ds 

1 

J' 

which is our approximate single-field equation for the FF. 



-|- {h.-o. linear and nonlinear terms} + O ( — ) =0, (61) 
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